Ground-state phase diagram of the spin-1/2 square-lattice 31-32 model with plaquette 

structure 
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Using the coupled cluster method for high orders of approximation and Lanczos exact diagonal- 
ization we study tire ground-state phase diagram of a quantum spin-1/2 J1-J2 model on the square 
lattice with plaquette structure. We consider antiferromagnetic (Ji > 0) as well as ferromagnetic 
(Ji < 0) nearest-neighbor interactions together with frustrating antiferromagnetic next-nearest- 
neighbor interaction J2 > 0. The strength of inter-plaquette interaction A varies between A = 1 
(that corresponds to the uniform J1-J2 model) and A = (that corresponds to isolated frustrated 
4-spin plaquettes). While on the classical level (s — 00) both versions of models (i.e., with ferro- 
and antiferromagnetic Ji) exhibit the same ground-state behavior, the ground-state phase diagram 
differs basically for the quantum case s = 1/2. For the antiferromagnetic case (Ji > 0) Neel 
antiferromagnetic long-range order at small J2/J1 and A > 0.47 as well as coUinear striped an- 
tiferromagnetic long-range order at large J2/J1 and A > 0.30 appear which correspond to their 
classical counterparts. Both semi-classical magnetic phases are separated by a nonmagnetic quan- 
tum paramagnetic phase. The parameter region, where this nonmagnetic phase exists, increases 
with decreasing of A. For the ferromagnetic case (Ji < 0) we have the trivial ferromagnetic ground 
state at small J2/IJ1I. By increasing of J2 this classical phase gives way for a semi-classical pla- 
quette phase, where the plaquette block spins of length s = 2 are antiferromagnetically long-range 
ordered. Further increasing of J2 then yields coUinear striped antiferromagnetic long-range order 
for A > 0.38, but a nonmagnetic quantum paramagnetic phase A < 0.38. 



I. INTRODUCTION 

The spin-1/2 quantum Heisenberg antiferromagnet 
with nearest-neighbor (NN), Ji > 0, and next-nearest- 
neighbor (NNN) bonds, J2 > 0, on the square lattice has 
attracted much interest (see, e.g., Refs.lll- [25l) as a canon- 
ical model to study the interplay between frustration 
and quantum fluctuations. In particular, the quantum 
phase transitions inherent in this model as well as the 
nature of its quantum paramagnetic phase in the region 
0.4 < J2/J1 ^ 0.6 is a matter of intensive debate. The 
Neel antiferromagnetic (NAF) long-range order (LRO) at 
small J2/J1 and the coUinear striped antiferromagnetic 
(CAF) LRO at large J2/J1 correspond to their classi- 
cal counterparts, however, the sublattice magnetization 
is reduced by quantum fluctuations. 

Motivated by recent investigations on quasi-two- 
dimensional frustrated magnetic materials with ferro- 
magnetic (FM) NN bonds e.g., Pb2VO(P04)2r^ — 



SrZnVO(P04)2 



29,31-33 



BaCdVO(P04)2r^'^'^ the FM 



Ji- J2 model ( Ji < 0, J2 > 0) has recently been 
studiedii^i^"— Although on the classical level the ground 
state (GS) phase diagrams of both variants of the Ji- J2 
model are quite similar, the GS properties of the quan- 
tum model with Ji < are basically different from those 
for Ji > 0. Currently, the GS properties of the FM model 
around J2/I Ji| = 0.5 are under controversial debate. On 
the one hand in Refs. [s^ and 13 arguments for a nematic 
phase separating the classical FM and the semi-classical 
CAF phases are given, on the other hand in Refs. l43l and 
I45I an intermediate quantum paramagnetic GS phase is 
not found (or it exists in a very small parameter region 



around J2 ^ 0.4|Ji|, only). Hence, a direct first-order 
transition between the FM GS and the GS with CAF 
LRO at J2 ^ 0.4| Ji| could take place4^ 

Quantum phase transitions can also occur by competi- 
tion between antiferromagnetic (AFM) NN bonds of dif- 
ferent strength, i.e., without frustration. One example 
is the local singlet formation in dimerized Heisenberg 
models, see e.g. Ref. 13 and references therein. As a 
certain extension to dimerized models the square-lattice 
Heisenberg model with a plaquette structure has been 
considered f^i^^— where local quadrumer singlet forma- 
tion can destroy NAF LRO. 

However, for the description of real materials a modifi- 
cation of the J1-J2 model might be necessary. For exam- 
ple, for PbZnVO(P04)2 a spatially anisotropy model has 
been derived;^ whereas for (CuCl)LaNb207i^ a J1-J2 
model with additional plaquette structure was proposed 
in Ref. IH- This plaquette model, however, has been 
questioned recently by Rosner and coworkersi^ 

In this paper we consider the J1-J2 model on the square 
lattice with plaquette structure (see Fig. [1]) merging 
this way the two mechanisms to destroy magnetic LRO, 
namely frustration and local singlet formation. More- 
over, the model we will consider corresponds to that pro- 
posed in Ref. 13 for (CuC^LaNbsOr. The Hamihonian 
of our model is given by 

H = ^((5oa+<5laA) Ji ^ SjSj + J2 ^ SiSj , (1) 

" = V {{ij))a ) 

where the sum is taken over the nearest ((•••)) and next- 
nearest (((■■■))) neighbors, and 5ab is the Kronecker sym- 
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bol. Wc consider FM and AFM NN bonds Ji = ±1 and 
AFM (i.e. frustrating) NNN bonds J2 > 0. The intra- 
(a ~ 0) and intcr-plaquette interaction (a = 1) differ by 
the factor A. The case A = 1 corresponds to the stan- 
dard uniform J1-J2 model on the square lattice, whereas 
the limit A = describes unconnected 4-spin J1-J2 pla- 
quettes. For the parameter A we consider the interval 
A G [0,1] this way interpolating between a system of iso- 
lated plaquettes and the uniform J1-J2 model. 




FIG. 1: The two-dimensional spin 1/2 J1-J2 finite square 
lattice of A'' = 32 sites (periodic boundary conditions) with 
plaquette structure (see Eq. ([T])); sohd lines - NN bonds Ji 
(thick) and AJi (thin), dashed lines NNN bonds J2 (thick) 
and AJ2 (thin). 

The model ([T|) has been studied previously by series 
expansion^ for AFM Ji > 0, and by bond-operator mean- 
field theory as well as a second-order perturbation theory 
in A^ for (FM) Ji < 0. Bearing in mind the intensive 
work on the standard J1-J2 model, see Refs. [1I-I25I and 
references therein, it seems to be desirable to discuss the 
much less studied plaquette J1-J2 model by methods be- 
ing alternative to those used in Refs. and this way 
to confirm or to question the findings of those papers. 

In this paper we consider both, the AFM ( Ji > 0) and 
the FM ( Ji < 0) version of the model ([1]), and we derive 
the GS quantum phase diagram for both versions of the 
model. We use the Lanczos exact diagonalization, see 
Sec. nil Al and the coupled cluster method (CCM), see 
Sec, nil B1 to analyze the GS of the model. Both methods 
are powerful general many-body methods and, particu- 
larly, the combination of both can lead to a consistent 
description of various aspects of quantum spin systems. 
We note that another important method for quantum 
spin systems, the quantum Monte Carlo method (QMC), 
cannot be used for the considered frustrated model, since 
it suffers from the minus sign problem. 



II. CLASSICAL GROUND-STATE PHASES 

It is well-known that the GS for the classical case 
(s — ?> 00) of the uniform J1-J2 model (i.e., A = 1) has 
three phases: The NAF phase for J2 < Ji/2 and Ji > 0, 
the FM phase for J2 < — Ji/2 and Ji < 0, and a phase 
consisting of two interpenetrating Ncel-ordercd square 
lattices. The angle between the directions of the two 
interpenetrating Neel states is arbitrary in the classi- 
cal limit, whereas due to the order- from- disorder effect 
coUinear states are favored by fiuctuations, i.e., the CAF 
is realized. For A < 1 to the best of our knowledge the 
classical phase diagram was not discussed so far. By 
contrast to the uniform model, two additional plaquette 
phases appear for A < 1, see Fig.[2l namely the plaquette 
antiferromagnetic (PAF) phase separating the NAF from 
the CAF phase for Ji > and the plaquette ferromag- 
netic (PFM) between the FM and the CAF phases for 
Ji < 0. These plaquette phases have some relation to the 
CAF phase: They also can be understood as a system of 
two interpenetrating Neel-ordered square lattices, where 
again the angle between the directions of the two inter- 
penetrating Neel states is arbitrary in the classical limit. 
However, as an elementary unit (a 'lattice site') of the 
two interpenetrating square lattices does not act a single 
spin (site), but rather a 4-spin-(4-site-)plaquette carrying 
strong bonds. And again, due to the order- from- disorder 
effect coUinear states are favored by fluctuations, i.e., the 
collinear PAF or PFM phases are realized, respectively. 
A graphical illustration of these classical phases is given 
in Fig. [21 The energy of the classical plaquette phases 
reads iJpAF.PFM = ~ [ii-^i ~ 1-^2] — |AJ2, where the 
upper (lower) sign corresponds to PAF (PFM). Obvi- 
ously, the inter-plaquette interaction is due to A J2 only. 
The transition line between the CAF and the PAF (PFM) 
phases is given by J2 = Ji/(1 -f A) (J2 = — Ji/(1 -|- A)), 
the transition between the NAF and the PAF phases as 
well as the FM and the PFM phases are horizontal lines 
J2 = >/i/2 and J2 = — Ji/2, respectively. Except the 
trivial degeneracies, the CAF, the PAF and the PFM 
states are two-fold degenerate, since the stripes of paral- 
lel spins (CAF) or parallel plaquettes (PAF, PFM) can 
be arranged cither along the horizontal or the vertical 
direction. 



III. METHODS 

A. Exact diagonalization of finite lattices 

The Lanczos exact diagonalization method was suc- 
cessfully used to discuss the GS phases of the uniform 
s = 1/2 J1-J2 model (i.e., A = 1) using finite lattices of 
N = 16, 20, 32, 36 and 40 sites,^"^i^i^ However, the new 
classical phases, PAF and PFM, appearing due to the 
plaquette structure do not fit to the periodic boundary 
conditions of the finite lattices of TV = 20, 36 and 40 sites, 
and. therefore, we do not have the possibility to perform 
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FIG. 2: Left panel: The classical phase diagram of the two-dimensional J1-J2 square-lattice Heisenberg model with plaquette 
structure, see Eq. JT]). For the AFM model (Ji = -1-1) the GS phases CAF, NAF, and PAF exist. For the FM model 
(Ji = —1) the CAF GS phase also exists, but the NAF and the PAF phases are replaced by the FM and the PFM phases, 
respectively. Right panel: Illustration of the classical phases (CAF: coUinear striped antiferromagnetic order, PAF: plaquette 
antiferromagnetic order, NAF: Neel antiferromagnetic order, FM: ferromagnetic order, PFM: plaquette ferromagnetic order). 



a finite-size extrapolation used previously for the uni- 
form niodeli^i^i^ Using Jorg Schulenburg's spinpacl^ 
we have calculated the GS focusing on the finite lattice 
of = 32 sites shown in Fig. [T] for a large set of A and 
J2 values fixing Ji either to +1 or to —1. 

To analyze the GS magnetic ordering we have calcu- 
lated the spin-spin correlation function (spSi) as well as 
a general finite-size order paramete r^^i^^ 



1 



N-4 



{N - 4)2 



(2) 



that adds the total strength of the spin-spin correlation 
functions. Bearing in mind the strong intra-plaquette 
correlation functions for A < 1 we have excluded these 
correlations in the sum, this way increasing the weight of 
distant correlation functions. 



B. Coupled Cluster Method 

We first mention that the coupled cluster method 
(CCM) yields results directly in the thermodynamic limit 
N ^ 00. The CCM has been previously reviewed in sev- 
eral articles (e.g., for the AFM^ii^^il and for the FM 
J1-J2 model^i^ as well as for dimerized models^"—) 
and will not be repeated here in detail. For more general 
information on the methodology of the CCM, see, e.g., 
Refs.MM- 

The CCM is a quantum many-body method and is de- 
fined by a reference (or model) state |$) and a complete 



set of mutually commuting many-body creation opera- 
tors C/. For our model we choose the classical GSs, see 
Fig. [2I as reference states. It is convenient to perform an 
appropriate rotation of the local axis of the spins such 
that in the rotated coordinate frame the reference state is 
a tensor product of spin down states |$) = .... 
The creation operators are then the multispin creation 
operators = , ^t^j^^ ^t^t^t'''' where the in- 
dices k, . . . denote arbitrary lattice sites. 

The CCM parameterizations of the ket- and bra- GSs 
are given by (with Cy = (C/)+) 



(*| = ($|^e-^, s = i + J2SiCy. 

1^0 



(3) 



Using ($|C+ = = V/ ^ 0, C+ = 1, the 

orthonormality condition {^\CiC^\^) = 6ij, and the 
completeness relation ^ C/|(f>)((f>|C/ = 1 = |(f>)((f>| -f- 

^C/|$)($|C/ we get a set of non-linear and linear 

equations for the correlation coefficients Si and Si, re- 
spectively. The order parameter (sublattice magnetiza- 
tion) in the rotated coordinate frame is given by 



1 ^ 



(4) 



In the CCM the only approximation is the truncation of 
the expansion of the correlation operator S and S. We 
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use the well established LSUBm scheme, where all multi- 
spin correlations on the lattice with m or fewer contigu- 
ous sites are taken into account. The number of these 
configurations can be reduced using lattice symmetry 
and conservation laws, but it is increasing very rapidly 
with m. In the highest order of approximation considered 
here, LSUBIO, we have, for example, 180957 configura- 
tions for the coUincar stripe reference state (CAF) and 
219446 for the antiferromagnetic plaquette state (PAF), 
i.e., finally 180957 or 219446 corresponding coupled non- 
linear equations which have to be solved numerically. 



Since the LSUBm approximation scheme becomes ex- 
act for m — ^ oo, we can improve our results by extrapolat- 
ing the "raw" LSUBm data to m — >■ oo. There is ample 
empirical experience regarding how one should extrap- 
olate the magnetic order parameter M{m) for systems 
with quantum phase transition between magnetically or- 
dered and disordered GS phases. Following Refs. [iBl - 
we use M(m) = bo + &i(l/m)i/2 + h^iXlmfl'^ to 
extrapolate to m — > cxd. As a rulei^"— the lowest level 
of approximation, LSUB2, is excluded from extrapola- 
tion. For the particular 4-spin plaquette model consid- 
ered here, we expect LSUB2 to be especially poor, be- 
cause its cluster size (2 spins) is smaller than the unit cell 
(4 spins) of the model. Hence, we use LSUB4, LSUB6, 
LSUB8 and LSUBIO data for the extrapolations. 
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FIG. 4: ED data for the finite -size order parameter vtl' jrr^p^i 
as defined in Eq. ((2} versus J2 (AFM Ji = for various 
values of A for the finite lattice of A'^ = 32 sites. For conve- 
nience rr? is scaled by its value for the ferromagnetic state 
rr^FM = 0.25. 



IV. RESULTS FOR THE QUANTUM s = 1/2 
MODEL 



A. Antiferromagnetic nearest-neighbor exchange ,/i 
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FIG. 3: ED data for the spin-spin correlation functions (soSi) 
versus J2 (AFM Ji = H-1) for two values of A for the finite 
lattice of A*' = 32 sites. Except the site indeces and i cor- 
responding to Fig. [T]we also give the square of separation B? 
of the sites and i. Note that (sqSi) and (S0S4) are intra- 
plaquette correlation functions, whereas the other ones are 
inter-plaquette correlation functions. 



Wc first consider the AFM case, set Ji — -f 1, and start 
with the discussion of ED data, see Figs. |3] and ID Cor- 
responding to the three classical phases there are three 
regimes in the quantum model. The spin-spin correla- 
tion functions (sgSi) presented in Fig. [3] for A = 0.5 and 
A = 0.7 arc quite strong for small and large J 2 thus in- 
dicating semi-classical magnetic LRO, where the signs of 
(soSi) fit to the classical phases NAF and CAF. In an in- 
termediate regime, around J2 = 0.5, the inter-plaquette 
spin-spin correlations are weak. These three different 
regimes are also well seen in Fig. 21 where the finite-size 
order parameter vn?, see Eq. ([5]), is shown. In particular, 
the m? data yield clear evidence for regions with weak 
magnetic order. The widths of those regions increase 
drastically with decreasing A. Bearing in mind that for 
A = 1 the region of rnagnetic disorder is 0.4 < J2 < 0.6, 
see, e.g., Refs. [Uli^llilaillli, andH our ED data for 
A^ = 32 suggest that there is no semi-classical antifer- 
romagnetic LRO of plaquette type (PAF). Moreover, for 
small values of A < 0.5 the finite-size order parameter 
is small in the whole range of J2 values indicating the 
absence of any magnetic LRO. 

Now we discuss the CCM results {N 00). Fig. [5] 
provides the CCM GS energy per site eo for A = 0.5 
and A = 0.7 compared with the corresponding ED data. 
Obviously, the CCM and the ED agree well. To discuss 
magnetic LRO we consider the CCM-LSUBm magnetic 
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order parameter M (sublattice magnetization) extrapo- 
lated to m — > oo, see Eq. that is depicted in Fig. [HI 
To illustrate the quality of the used extrapolation for the 
order parameter M of the 'raw' LSUBjti data to the limit 
m — >■ oo we also show, as an example, corresponding plots 
for the A = 0.7 in Fig. [71 It is obvious that the LSUBm 
data are well fitted by the applied extrapolation function. 
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FIG. 7: Illustration of the extrapolation (solid lines) of the 
CCM-LSUBm data (symbols) for the magnetic order pa- 
rameter M for A = 0.7, Ji = -fl and specific values of 
J2 corresponding to the NAF, PAF, CAF reference states. 
For the extrapolations to m ^ 00 according to M{m) = 
bo + 6i(l/m)i/2 -I- 62(1/^)^/2 we have used LSUBm data for 
771 = 4,6,8, 10. 



FIG. 5: COM LSUBIO data for the GS energy eo (AFM 
Ji = +1) for A = 0.5 and A = 0.7 compared with corre- 
sponding ED data (A*' = 32). The CCM results correspond to 
the NAF reference state (small J2), the PAF reference state 
(intermediate J2), and the CAF reference state (large J2). 
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FIG. 6: Extrapolated CCM data for the magnetic order 
parameter (sublattice magnetization) A'l for the AFM case 
(Ji = -1-1) and for various values of A. The CCM reference 
state is the NAF state (left hand side, small J2) and the CAF 
state (right hand side, large J2). The extrapolated CCM data 
are obtained using LSUBm with m = 4, 6, 8, 10 and the ex- 
trapolation rule M{m) ^ bo + 6i(l/m)^/2 + b2{l/mf/^. 

First we notice that the extrapolated order parameter 
using the PAF reference state for all relevant parameter 
sets of A and J2 vanishes (see also Fig.[71), i.e., there is no 



PAF magnetic LRO in the quantum model. This find- 
ing is in agreement with the ED data for discussed 
above (and see again Fig. |4]) . The range of stability of the 
semi-classical NAF and CAF phases is visible in Fig. [51 
We define the quantum critical points J2^(A) and J^^iX) 
as those points, where the extrapolated CCM magnetic 
order parameter M vanishes. For A = 1 (uniform model) 
we find J|i = 0.447 Ji and Jf = 0.586Ji which is in 
agreement with previous CCM predictions J^"— As ex- 
pected, with decreasing of A the region of J2 values with- 
out magnetic LRO increases, see Fig. El Collecting the 
data for the quantum critical points and for var- 
ious A we get the GS phase diagram of the quantum 
model as shown Fig. |8l Our phase diagram is in good 
agreement with the only available corresponding one of 
Ref. i. Hence, one can argue that the phase diagram is 
basically correct. 

Concerning the order of the quantum phase transitions 
at and our data are in favor of a continuous 
transition at and a first-order transition at J^^ as it is 
was discussed for the uniform modeli^'^iii'^i The scenario 
of a first-order transition between the QPM and the CAF 
phases is supported by (i) the kink-like behavior in eo, 
see Fig. |5l (ii) the steep fall in the CCM order parameter 
M near J2^, see Fig. |S1 and (iii) the jump- like behavior 
of the spin-spin correlations functions, see Fig. |3l and of 
the finite-size order parameter, see Fig. [H 

For the unfrustrated case J2 = we can compare 
our result Ac = 0.47 (and see Fig. [5]) for the crit- 
ical A where the NAF LRO breaks down with sev- 
eral previous results, namely Ac = 0.555 (Ising series 
expansion^) Ac = 0.6 (exact diagonalization^) Ac = 
0.5485 (QMGi^i^), Ac = 0.5491 ... 0.5513 (contractor 
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FIG. 8: The COM GS phase diagram of the model O for the 
AFM case (Ji = +!)■ NAP denotes the semi-classical phase 
with Neel LRO, CAP the semi-classical phase with coUinear 
striped antiferromagnetic LRO, and QPM the magnetically 
disordered quantum paramagnetic phase. Por comparison we 
also show the classical transition lines (thin dashed lines), cf. 
Fig-H 



renormalization expansion^), and Ac = 0.4822 (real 
space renormalization group approach^). Likely, the 
QMC result is most accurate. Our result is in reasonable 
agreement with that result, but slightly overestimates the 
stability region of NAF LRO. 

Let us briefly discuss another limit of the model, 
namely the limit of large J2- As discussed above in that 
limit the system splits into two interpenetrating square- 
lattice Heisenberg antiferromagnets (where J2 plays the 
role of the AFM NN bond), which are Neel-ordered for 
A = 1. Hovewer, for A < 1 each of the two interpenetrat- 
ing square lattices carries now a staggered arrangement of 
J2 and A J2 bonds, which corresponds precisely to the so- 
called J- J' model discussed in Refs. [EflillleiHi-IlO. The 
QMC estimate of the critical A is Ac = 0.397, see Ref . Isol. 
For J2 = 1 we find from our CCM data a critical value of 
Ac = 0.301 (see Fig. |8]). Increasing of J2 beyond J2 = 1 
(not shown in Fig. HI) yields a steep increase of the criti- 
cal line as already indicated by the last two data points 
near J2 = 1 in Fig. [8] In the limit of large J2 the critical 
Ac is even smaller as for J2 = 1. For example at J2 — 5 
we obtain Ac — 0.375 that is in good agreement with the 
QMC result for the critical Ac of the J- J' model. Note 
that our value of Ac deviates from an early CCM result 
Ac = 0.316.— The difference is related to the fact that 
we use here (i) a higher approximation level (LSUBIO) 
and (ii) an improved extrapolation in comparison with 
Ref. leO- 



B. Ferromagnetic nearest-neighbor exchange Ji 

We now consider the FM case and set Ji = —1. Al- 
though, the classical phase diagrams for Ji = — 1 and 



FIG. 9: ED data for the spin-spin correlation functions (soSi) 
versus J2 (PM Ji — —1) for two values of A for the finite 
lattice of A'' = 32 sites. Except the site indeces and i cor- 
responding to Fig. [1] we also give the square of separation 
of the sites and i. Note that (sosi) and (S0S4) are intra- 
plaquette correlation functions, whereas the other ones are 
inter-plaquette correlation functions. 




FIG. 10: ED data for the finite-size order parameter m /rupj^i 
as defined in Eq. ([2]) versus J2 (FM Ji = —1) for various 
values of A for the finite lattice of A'^ = 32 sites. For conve- 
nience is scaled by its value for the ferromagnetic state 
= 0.25. 



Ji = +1 arc identical, see Fig. [21 we know from the 
uniform model (i.e., A = 1) that the phase diagrams in 
the quantum case s = 1/2 are basically different J^i^"— 
In particular, it was found for Ji = — 1 that the region 
of semi-classical CAF LRO extends up to much smaller 
values of J2 compared with the case of Ji = 
Hence, we expect also for A < 1 basic differences between 
Ji = +1 and Ji = —1. 

We start again with ED results for the spin-spin corre- 
lation functions (soSi), see Fig. [HI and the finite-size order 
parameter defined in Eq. Q , see Fig. (TUl The three 
classical phases lead again to three different regimes in 
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FIG. 11: CCM-LSUBIO (thick lines) and LSUB6 data (thin 
lines) for the GS energy eo (Ji — —1) for A = 0.5 and A = 0.7 
compared with corresponding ED data (A^ = 32). The CCM 
results correspond to the PFM reference state (intermediate 
J2) and the CAF reference state (large J2). The black solid 
line is the exact FM GS energy. 



the quantum model. The trivial FM state at small J2 is 
of course also present in the quantum model. It gives way 
at J|^(A) for an intermediate singlet state. This transi- 
tion at is clearly seen by jumps in (soSj) and m^. 
is smaller than the classical value J2 cZas = 0.5, and it 
depends on A for the quantum model. Moreover, there is 
a second jump-like behavior of the spin-spin correlation 
functions indicating the transition from the intermediate 
regime to the CAF regime. We can use the positions of 
these jumps to extract the transition points and 
between the regimes from our ED calculations, see the 
discussion of Fig. [1^ given below. 

The main difference in comparison to the AFM model 
(Ji = -1-1) concerns the intermediate regime. From Fig. [5] 
it is evident that the inter-plaquette correlation functions 
(sqSio) and (S0S22), see Fig. [U are not small. Hence, 
the intermediate phase might be long-range ordered for 
FM Ji. Indeed, in the finite-size order parameter 
shown in Fig. [TU] an intermediate regime is clearly vis- 
ible for A < 0.5, where the finite-size order parameter 
m? in this regime is much larger than that for Ji = +1 
(cf. Fig. 21) . The possible appearance of magnetic LRO 
can be understood on the basis of the classical PFM 
state (see also Sec. |ll|: The elementary unit is the 4- 
spin-plaquette (with strong bonds). Due to the strong 
FM intra-plaquette bond Ji the plaquette carries an ef- 
fective block (composite) spin s = 2. These block spins 
interact via A J2 and form effectively two interpenetrating 
square-lattice s = 2 Heisenberg antiferromagnets. 

Now we pass to the CCM results {N ^ 00). Fig. [TT] 
provides the CCM GS energy per site eo for A = 0.5 and 
A ~ 0.7 in LSUB6 and LSUBIO approximation compared 
with the corresponding ED data. The agreement between 
the CCM and the ED is good. The three regimes, FM, 



FIG. 12: GS phase diagram of the plaquette model of Eq. ([T]) 
for the FM case (Ji = —1). The transition lines between 
the FM and the PFM phases as well as the PFM and the 
CAF/QPM phases are obtained by (i) ED by using the jumps 
in the ED data (black lines), see Figs. and 1 101 and (ii) by 
CCM from the intersection points between the GS energies 
using various reference states for the CCM calculation (blue 
lines), see text. The transition lines between the CAF and the 
QPM phase (red line) correspond to those parameter values 
(A, J2), where the extrapolated CCM order parameter van- 
ishes (and see Fig. I14p . For comparison we also show the 
classical transition lines (magenta line), cf. Fig. [2] 
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FIG. 13: CCM magnetic order parameter M for the FM 
case (Ji = —1) versus J2 using the PFM and the CAF 
states as CCM reference states for A = 0.5 and A = 0.7. 
The extrapolated CCM results are obtained using LSUBm 
with m — 4, 6, 8, 10 and the extrapolation rule M{m) = 
bo + bi{l/my/^ + b2{l/mf/^ . 



PFM, and CAF, are clearly seen in the energy data. For 
the CCM estimate for the transition point between 
the FM state and the PFM regime we use the intersection 
points between the FM energy and the CCM-LSUBIO as 
well as the LSUB6 energy calculated with the PFM ref- 
erence state. Both results for almost coincide, see 
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FIG. 14: Magnetic order parameter M for the FM case 
(Ji = —1) versus A for J2 = 1.0 (filled triangles), 0.7 (circles) 
and 0.5 (lines) using the PFM and the CAF states as CCM 
reference states. The extrapolated CCM results are obtained 
using LSUBm with m = 4, 6, 8, 10 and the extrapolation rule 
M(m) = bo + &i(l/m)^/2 + b2(l/m)3/^ 



Fig. [T^l Moreover, we find an excellent agreement be- 
tween the ED and CCM results for J|^. Concerning the 
order of the phase transition at we have clear evidence 
for a first-order transition as indicated by (i) the jump- 
like behavior of the spin-spin correlations functions, see 
Fig. m and of the finite-size order parameter, see Fig. [101 
(ii) the kink-like behavior in eo, see Fig.[TTJ and (iii) the 
jump in the CCM order parameter M at J|^, see Fig.[T3l 

To determine the transition point J!f' between the 
PFM and the CAF regimes we are faced with the prob- 
lem that near this transition no CCM-LSUBIO solutions 
for the CAF reference state could be found. This ob- 
servation, that for high orders of CCM approximation in 
the vicinity of an intersection point of two CCM energy 
curves belonging to two different GS phases no solution 
of the large set of coupled nonlinear ket equations can be 
obtained, is often found when dealing with strongly frus- 
trated systems, sec, e.g., Refs. [isl andlisl. However, from 
Fig. E] it is obvious that the parameter range, where so- 
lutions for lower levels m of LSUBm approximations can 
be found, is much larger than for LSUBIO. Moreover, 
the LSUBm data for cq for various m are very close to 
each other. Hence, we can use the intersection points 
for lower levels of CCM-LSUBm approximation to de- 
termine the second transition point J|^. Thus, wc find 
intersection points for LSUB8, LSUB6, and LSUB4 for 
A > 0.7, A > 0.5, and A > 0.25, respectively. In addition, 
we can also take benefit from the almost straight behav- 
ior of eo(J2) curves, see Fig. [Til which allows a reliable 
extrapolation of eo( J2) until a hypothetical intersection 
point, cf. also Refs. and H^. This gives finally vari- 
ous sets of J|^(A) data, which, however, agree well with 
each other, cf. Fig. [T^l Only for A < 0.3 a slight differ- 
ence is visible. Nevertheless, it is necessary to mention 



that the CCM estimate for J|^(A) becomes less reliable 
for smaller values of A due to the increasing distance be- 
tween the hypothetical intersection points and the last 
data points for which LSUBm solutions for the CAF ref- 
erence state could be found. The results for (A)and 



0.5 




FIG. 15: Illustration of the extrapolation (solid lines) of 
the CCM-LSUBm data (symbols) for the magnetic order 
parameter M for A = 0.5, Ji = —1 and specific values 
of J2 corresponding to the CAF and PFM reference states. 
For the extrapolations to m — ^ cx) according to M{m) = 
ho + 6i(l/m)^/^ + 62(l/m)^''^ we have used LSUBm data for 
m = 4,6,8, 10. 

J|^(A) obtained by ED and CCM are collected in the 
phase diagram presented in Fig. 1121 It is obvious that 
both, ED and CCM, yield very similar values for and 
^2 ■ 

The question about magnetic LRO in the various 
regimes wc address next by analyzing the CCM mag- 
netic order parameter M for the PFM and CAF regimes. 
(Remember that in the trivial FM phase we have al- 
ways M = 1/2.) In Fig. [13] we show M versus J2 for 
A = 0.5 and A = 0.7. Obviously, the order parameter 
in the PFM regime is non-zero and even larger than in 
the CAF regime. This observation is in agreement with 
our ED results, cf. Fig. [TOl and can be related to the 
effective block-spin s = 2 model discussed above. The 
order parameter in the CAF is only weakly dependent 
on J2, but its magnitude depends on A. As discussed 
above we do not get CAF CCM solutions for LSUBm 
for higher values of m near the transition to the PFM 
phase. Nevertheless, the curves presented in Fig. [T31 in- 
dicate that there is likely a direct first-order transition 
at between semi-classical phases with PFM and CAF 
magnetic LRO for A > 0.4. 

Next we fix J2 and consider the behavior of the order 
parameter M in dependence on A, see Fig. [TH For val- 
ues of J2 ^ 0.55 the PFM regime is not relevant. As 
explained in Sec. IIV Al in the limit of large J2 our pla- 
quette model corresponds to the J-J' model, where at 
A = Ac = 0.397 (QMC result, see Ref.llO) a second-order 
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transition to a QPM phase takes place. This behavior is 
clearly seen in our CCM data for M(A) for J2 = 1-0 and 
J2 = 0.7 shown in Fig. [H We find that for J2 > 0.6 the 
extrapolated CCM order parameter M vanishes contin- 
uously at a critical value A = Ac(J2) defining a second- 
order transition between a semi-classical phase with CAF 
magnetic LRO and the QPM phase, that is depicted as a 
red line in Fig. [T^ Interestingly, Ac practically docs not 
depend on J2 and we have Ac ~ 0.383 for 0.6 < J2 < 1- 
This value is close to the QMC valued Ac = 0.397 valid 
for J2 — >■ cx) and also to Ac = 0.375 obtained for the AFM 
model (Ji = -1-1) for J2 = 5, see Sec. lIVAl 

For J2 = 0.5 the situation is different, since both, the 
PFM and CAF regimes, play a role. For A < 0.494 we are 
inside the PFM regime. Obviously, the PFM order pa- 
rameter is large and the variation with A is weak. Within 
the CAF regime we are again faced with the problem that 
the CAF LSUB7T1 solutions for higher m terminate be- 
fore meeting the corresponding PFM LSUBm solutions. 
Hence, the critical line Ac( J2) (red vertical line in Fig. fT^ 
terminates before meeting the transition line ^2^ (A), and 
we cannot give a secure statement on the continuation of 
the critical line Ac(J2) towards lower values of J2- 

Finally, we illustrate the quality of the used extrapola- 
tion for the order parameter M of the 'raw' LSUBm data 
to the limit m — 5- 00 for A = 0.7 in Fig. [131 It is obvi- 
ous that the LSUBto data are well fitted by the applied 
extrapolation function. 

Let us briefly discuss the relation of our results to the 
schematic phase diagram previously presented in Ref. |48| . 
As stated in Ref. |4^ the approximations used there 
(bond-operator mean-field theory as well as a second- 
order perturbation theory in A) may be not reliable for 
large A ^ 1 and large J2 ~ 1. The main features of our 
phase diagram agree well with those presented in Ref. IH. 
For A — > we obtain the same transition point J2^ = 0.5 
between the QPM and the PFM phases. However, for 
A 1 and J2 ~ 1 some differences appear. For instance, 
we do not find a nematic phase for A < 1, as it was dis- 
cussed in Ref. Note that the absence of the nematic 



phase is in agreement with the recent findings for the 
uniform model, i.e. at A = iM^^ Moreover, the second- 
order transition between the QPM and CAF phases for 
large J2 is found in Ref. [H at A ~ 0.5, but it should be 
at A ?» 0.4. 



V. SUMMARY 

Inspired by a recent investigation of a frustrated two- 
dimensional Heisenberg model proposed to describe the 
magnetic properties of (CuCl)LaNb207, see Refs.ls^and 
m, we investigate the GS phase diagram of the spin- 
1/2 J1-J2 Heisenberg model on the square lattice with 
plaquette structure. The 4-site plaquettes carrying the 
(strong) intra-plaquette bonds Ji and J2 are coupled to 
each other by (weaker) inter-plaquette bonds AJi and 
^•h: < A < 1. The parameter A can also be thought of 
modeling a distortion of the underlying square lattice. 

We consider AFM (Ji > 0) as well as FM nearest- 
neighbor exchange coupling (Ji < 0). Except the 
phases with magnetic LRO (FM, Neel, and collinear 
striped AFM) and the non-magnetic quantum param- 
agnetic phase known from the standard spin-1/2 J1-J2 
model we find for FM Ji a new plaquette phase showing 
antiferromagnetic long-range of s = 2 block spins associ- 
ated with 4-spin plaquettes. For the AFM model ( Ji > 0) 
the region of the quantum paramagnetic phase is signif- 
icantly increased for A < 1 compared to the standard 
model, thus increasing the prospects of finding a magnet- 
ically disordered low-temperature phase in real magnetic 
materials, where the exchange pattern may deviate from 
the standard J1-J2 model, see, e.g., Refs. IssIItiI and[72l. 
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